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This paper presents the implementation of a complex fractional order 
proportional integral derivative (CPID) and a real fractional order PID 
(RPID) controllers. The analysis and design of both controllers were carried 
out in a previous work done by the author, where the design specifications 
were classified into easy (case 1) and hard (case 2) design specifications. 
The main contribution of this paper is combining CRONE approximation 
and linear phase CRONE approximation to implement the CPID controller. 
The designed controllers-RPID and CPID-are implemented to control 
flowing water with low pressure circuit, which is a first order plus dead time 
system. Simulation results demonstrate that while the implemented RPID 
controller fails to stabilize the system in case 2, the implemented CPID 
controller stabilizes the system in both cases and achieves better transient 
response specifications. 
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1. INTRODUCTION 

Fractional calculus is an extension to ordinary calculus by extending the orders of differentiation 
and integration to noninteger numbers. This mathematical concept was utilized in system modelling and 
control. In control, several fractional controllers were designed for several types of systems [1], [2]. As a 
further extension, complex fractional calculus was introduced as an extension to the real fractional calculus, 
where the orders of differentiation and integration can be complex numbers [3]. In [4], definitions and 
theorems were presented for complex fractional calculus mathematically. Complex fractional calculus was 
utilized to introduce models that describe viscoelastic materials [5], to model drug resistance in human 
immunodeficiency virus (HIV) infection [6], [7], and to present a new mathematical model for the atrial 
fibrillation [8]. Cois et al. [9] proposed a tool to model and study state-space with complex order. A complex 
fractional calculus was also applied to solve certain mathematical problems, such as fractional boundary 
problems [10]. 

In control theory and applications, fractional calculus has been well utilized to design fractional 
order controllers. Since the conventional proportional integral derivative (PID) controller dominates other 
controllers [11], [12], the most common fractional order controller is the fractional order PID controller (also 
called PID! controller), proposed by I. Podlubny [13], [14]. It is a generalization of the conventional PID 
controller, where the integer order derivative and integral actions are replaced by fractional order derivative 
and integral actions. Some of the recent works that utilize the PID" controller are presented in [15]-[17]. 
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This generalization was also utilized in the differentiation and integration actions of the conventional fuzzy 
logic controller, resulting in a fractional order fuzzy logic controller [18], [19]. For nonlinear controllers, 
such as sliding mode controllers (SMC), fractional calculus was introduced in the sliding surface by taking 
the fractional order derivative and/or integral of the state variables [20]-[27]. Hybridization between fuzzy 
logic and fractional order sliding mode control can be achieved, where the resultant controller is called Fuzzy 
fractional order sliding mode controller (FFOSMC) [28]-{32]. All the previously mentioned works and other 
works in the literature lack the investigation of the effect of introducing complex fractional calculus into 
control action. In [33], the analysis of real fractional order PID (RPID) and complex fractional order PID 
(CPID) controllers were carried out, and both controllers were designed to control a first order plus dead time 
(FOPDT) system. In the design process, the specifications were classified into two cases, easy and hard. This 
classification is adopted to demonstrate the need for the CPID controller to overcome the limitations of the 
RPID controller. This paper implements the RPID and CPID controllers designed in [33]. The RPID 
controller is implemented using CRONE approximation-as most works in the literature; however, the 
contribution of this paper is the implementation of the CPID controller; it proposes combining CRONE 
approximation and linear phase CRONE approximation to obtain an acceptable approximation of this 
controller. 

The rest of the paper is organized: section 2 presents the problem statement; section 3 presents the 
model of the flowing water with low pressure circuit (FWLPC), which is the plant to be controlled; section 4 
presents the implementation of the RPID and CPID controllers using CRONE approximation and CRONE 
approximation/linear CRONE approximation, respectively; the results and their discussion is presented in 
section 4; and the conclusions is drawn in section 5. 


2. PROBLEM STATEMENT 

This paper addresses the utilization of complex fractional calculus in control theory. It is concerned 
with implementing the CPID controller that was designed in [33]. The suggested method is to combine 
CRONE approximation and linear phase CRONE approximation to approximate the real 
differentiator/integrator and imaginary differentiator/integrator factors of the controller, respectively. 


3. MATHEMATICAL MODEL OF THE FWLPC 
The LPRWC-the process to be controlled-is a FOPDT system; its transfer function is 


P(s) = —~e-#s (1) 


TS+1 


where L = 50s (time delay), k = 3.13 (gain), and T = 433.33s (time constant) [33], [34]. 


4. CONTROL SYSTEM IMPLEMENTATION 

The RPID and CPID controllers that are implemented in this paper were designed in [33], where the 
design specifications (desired gain crossover frequency w, and desired phase margin m) were classified into 
two cases; in case 1, the design specifications are easy, while in case 2, the design specifications were made 
hard by increasing both w, and Qm- 


4.1. Implementation of the RPID controller 
The RPID control law u(t) is 


u(t) = Kpe(t) + K,D~e(t) + KpD“e(t) (2) 
and the corresponding transfer function is 


Cr(s) = TS = Kp + Ki + Kps" => Kp + Kıs ~? + Kps” (3) 


where Kp, Kı, Kp E R and A and u E R* are the parameters of the controller were designed in [33] and are 


shown in Table 1. To implement Cg(s), it is needed to implement the fractional-order integrating action s72 


and fractional-order differentiating action s”. In this paper, CRONE approximation (sometimes called 
Oustaloup filter) is used. It is given by [35]. 
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Vis N St@z,k 

S K IR- POR (4a) 
2k-1-y 

Wz k = Wu) N (4b) 
2k-1+y 

Wp, k = W(@y,) N (4c) 

K =o (4d) 

o= [oe (4e) 


where [w;, wp] are the frequency range on which the approximation is valid. The frequency range is taken as 
[wn wn] = [0.01, 100] rad/s. Using the values of the RPID controller parameters given in Table 1 and 
(4a)-(4e) to approximate s72 and s” in (3), Cg (S) becomes 


690.258+2.012x10457+1.737X10°s©+3.978X10°5°+3.211X10°s*+7.689X104534+6925s2+2155+2.446 
58499.94574+1927s©+1.182X104s°+2.042 X10454+1.193 X1045341965s524+102.95+1.039 
(5) 


Ca(s) = 


for case 1, and 


Ca(s) ~ 
6.567X10*58+2.616x10°S7+1.979X1075°+4.785X107 5°+3.26X107 s*+7.816X10°s? +6.349X10°s2+3.0291045+545.7 
584254.1574+1.202 X10456+1.824X10°s°5+7.707X10°54+1.112X10553+4.47 X105s2+5.76X1045+1382 


(6) 
for case 2. 
Table 1. RPID and CPID controller parameters 
Parameters 
Case Controller 
Kp Kı Kp a B 0 $ 
1 RPID 0.09 0.02 5.22 0.88 0.31 
CPID 0.76 0.16 21.54 0.04 1.19 0.21 1.56 
2 RPID 17.40 0.07 29.30 1.76 0.64 
CPID 0.94 20.65 31.92 0.48 -1.42 0.85 0.73 
4.2. Implementation of the CPID controller 
The CPID control law u(t) is 
u(t) = Kpe(t) + KıD CHP e(t) + KpD? tt e(t) (7) 


and the corresponding transfer function is 


U(s) 1 ; 1 ; 1 
Cc(s) =o Kp + Ki -g + Kps? = C(s) = Kp + Ki zg + Kps? Y? = Kp + Kı cag t 


Kps? st? = Kp + K;s~%s J? + Kps? st? (8) 


where £, p €R are the extra parameters that introduces the imaginary integrator s~/ and imaginary 
differentiator s/?. The seven parameters of the CPID controller were designed in [33] and are shown in 
Table 1. The fractional integrator s~* and fractional differentiator s? are approximated using (4a)-(4e). 
However, to implement (8), it is needed to find a valid approximation of the imaginary integrator s~/* and 
imaginary differentiator s/?. The sinusoidal transfer function of s/” is 


SY = (jW) = jY x wF = ez x (cos(y In w) + j sin(y Inw)) (9) 


[goy] = e> (10) 
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From (11), it can be seen that the phase is a linear function of log,,, i.e., if the phase is plotted in a 
logarithmic axis, it is a straight line with slope equals to y In 10. This is a significant feature of the imaginary 
differentiator/integrator, since the real fractional differentiator/integrator has constant phase angle (i.e., its 
slope equals zero). In this paper, the imaginary differentiator/integrator are approximated by linear CRONE 
approximation, which relies on locating the transfer function poles and zeros of the approximating transfer 
function such as to give a linear phase curve. Unlike the CRONE approximation, in linear CRONE 
approximation, the spacing distance between the poles is not equal to the spacing distance between the zeros; 
this different spacing distance is responsible of giving the nonzero slope of the phase line. The linear CRONE 
approximation is [4], 


: 1t> 7 
s7 x~ TERN, TAE (12a) 
®py,k 
1 
wz k = We (b) "2 (12b) 
a k-N-+ 
pk = Wc(a) ? (12c) 


We = Jf @ Wn (12d) 


where [w, œn] are the frequency range on which the approximation is valid, a and b are the recursive 
coefficients of the poles and zeros, respectively. The relation between the imaginary order y and the recursive 
coefficients a and b is 


i) eee 
logiob =logioa (13) 


yln10 = 


Algorithm 1 shows a pseudo code to set the values of a, b, and N given the values of w);, wpn, and y. 
The frequency range of the linear CRONE approximation is also taken as [w,, wn] = [0.01, 100] rad/s (to 
give a fair comparison with the CRONE approximation of the RPID controller). Applying the algorithm 
given in Figure 1, the resultant values of a, b, and N are given in Table 2. 


Algorithm 1: A flow chart to set the values of a, b, and N given the values of w,, Wp, and y. 
1- Input @,, Wp, and y. 

2- Set a equals to a certain value >1 and set N=0. 

3- Calculate b from equation (13). 

4- N=N+1. 

5- Calculate w,, from (12b). 

6- Calculate W,, from (12c). 

7- If b<a and wk >W@p, go to 10. 

8- If b>a and Wpk > Wn, go to 10. 

9- Go to 4. 

10- Output the values of a, b, poles Wpkr ZELOS Wx, and N. 


Using the parameter values of the CPID controller in Table 1, (12a)-(12d), and Algorithm 1 to 
approximate sP and sJ® in equation (8), Cc(s), for case 1 and case 2, becomes 


Xt as! 


Co(s) = 27 pis! (14) 
for case 1, and 
Z2 cist 
Co2(s) ~ Sasi (15) 


for case 2, where a; and b; are given in Table 3 and c; and d; are given in Table 4. 
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Figure 1. Unit step response of the RPID and CPID control systems: case 1 


Table 2. Linear CRONE parameter values a, b and N and the resultant slope of phase line 


Case Imaginary differentiator/integrator a b N Slope of phase line 
1 gs J119 2 1.9280 4 0.2916 
g/156 2 2.1346 4 -0.4482 
2 sJC1-42) 2 1.5487 4 3.0505 
51073 2 1.8550 4 0.6358 


Table 3. Values of the coefficients a; and b; in (14) Table 4. Values of the coefficients c; and d; in (15) 


Coefficient Value Coefficient Value Coefficient Value Coefficient Value 
ao 43.1 bo 17.5 Co 1.8 do 0.9 
a 4664 bı 1776 Cy 194.1 dı 106.5 
ay 2.188 x 10° b 7.297 x 104 Cp 8113 dz 5079 
a3 5.8 x 10° bs 12.9 x 10° C3 1.8 x 105 dz 1.3 x 105 
a, 1x108 b, 2.3 x 107 C4 2.4 x 10° dy 2.1 x 10° 
as 1.1 x 10° bs 2.2 x 108 Cs 2.2 x 107 ds 2.2 x 107 
dg 9.8 x 10° be 1.5 x 10° C6 1.4 x 108 de 1.6 x 108 
ay 5.8 x 101° b; 7.1 x 10° Cy 6.7 x 108 d, 8.7 x 108 
ag 2.5 x 101 bg 2.4 x 101° Cg 2.3 x 10° dg 3.3 x 10° 
dg 8.1 x 1012 bo 6.3 x 101° Co 5.9 x 10° dy 9.4 x 10° 
aio 1.9 x 10% Dio 1.2 x 1014 Cio 1.1 x 101° dio 1.9 x 101° 
ait 3.4 x 10% bit 1.7 x 1014 Cu 1.7 x 107° diy 3x 10° 
a2 1.9 x 107 biz 1.8 x 101 ĝi? 1.9 x 107° dip 3.5 x 10° 
a13 4.3 x 10” bis 1.4 x 101 Cis 1.7 x 101° diz 3 x 101° 
a 3.16 x 10"? bas 8.4 x 101° Cia 1.1 x 10° di4 1.9 x 101° 
ais 1.71 x 101? Dis 3.7 x 101° Cis 5.7 x 10° dis 9.5 x 10° 
a16 6.8 x 1014 bie 1.2 x 107° C16 2.2 x 10° dig 3.4 x 10° 
a7 1.9 x 101 bız 2.8 x 10° C17 6.3 x 108 di7 9 x 108 
aig 4.2 x 107° big 4.9 x 108 Cis 1.3 x 108 dig 1.7 x 108 
aig 6.4 x 10° big 6x 107 Cio 2x107 dig 2.3 x 107 
a20 6.7 x 108 Doo 5.1 x 10° C20 2.2 x 10° doo 2.2 x 10° 
a21 4.7 x107 bot 2.8 x 10° C21 1.5 x 10° do, 1.3 x 10° 
a22 2.1 x 10° bo» 9316 C22 6943 dz2 5408 
a23 5.4 x 104 b23 160.4 C23 160.4 d23 114.6 
a24 582.8 boa 1 C24 1.4 dy, 114.6 


5. SIMULATION RESULTS AND DISCUSSIONS 

The implemented RPID and CPID controllers were used to control the FWLPC system given by (1), 
where the systems were simulated using MATLAB Simulink. The unit step responses for both systems were 
obtained for comparison purposes. Figure 1 shows these responses for case 1 and Figures 2-3 show these 
responses for case 2, and Table 5 shows the performance of both systems in terms of transient response 
specifications. The obtained results demonstrated that in case 1, the performance of the implemented CPID 
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controller is better than that of the implemented RPID controller by reducing both the rise time and steady 
state error, while in case 2 the RPID controller failed to stabilize the system while the CPID controller 
stabilized it. The significance of these results is that it insured the results that were obtained when the 
controllers were designed [33], where the RPID controller stabilized the system only in casel, while the 
CPID controller stabilized it in both cases and gave better design specifications. For both systems, the 
steady-state error is not zero; this is due to approximating the pure integrator in both the RPID and CPID 
controllers. 
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Figure 2. Unit step response of the RPID control system: case 2 


Liquid level (m) 


0 fi f fi f fi f fi f fi 
0 100 200 300 400 500 600 700 800 900 1000 


Time (s) 


Figure 3. Unit step response of the CPID control system: case 2 


Table 5. Transient response specifications of the RPID and CPID control systems 


RPID CPID 
Casel Case2 (Casel Case2 
Rise time (s) 70 unstable 62 90 
Percentage overshoot 7.5% 16% 
Steady-state error (m) 0.12 0.11 0.14 
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6. CONCLUSIONS 

In this paper, RPID and CPID controllers have been implemented to control a FOPDT system. The 
RPID controller has been approximated using CRONE approximation, and the CPID controller has been 
approximated by combining CRONE approximation and linear phase CRONE approximation. The following 
two conclusions are draw: i) adding extra implementation parameters has a significant impact on the order of 
the controller; 1i) the CPID controller outperforms the RPID controller from the stability and performance 
(transient response specifications) point of views. As a suggestion for future work, the steady-state error of 
both the RPID and CPID control systems can be investigated by identifying its cause and suggesting proper 
solutions to reduce or eliminate it. 
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